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A number of theories, spanning a wide range of mass scales, predict dark matter candidates that have lifetimes 
much longer than the age of the universe, yet may produce a significant flux of gamma rays in their decays 
today. We constrain such late decaying dark matter scenarios model-independently by utilizing gamma-ray line 
emission limits from the Galactic Center region obtained with the SPI spectrometer on INTEGRAL, and the 
determination of the isotropic diffuse photon background by SPI, COMPTEL and EGRET observations. We 
show that no more than ~5% of the unexplained MeV background can be produced by late dark matter decays 
either in the Galactic halo or cosmological sources. 

PACS numbers: 95.35.+d, 13.35.Hb, 12.60.-i 



I. INTRODUCTION 

Dark matter continues to live up to its name [1], despite 
accumulated evidence of its existence from observations of 
large-scale structure formation, galaxy cluster mass-to-light 
ratios, and galactic rotation curves. An attractive approach to- 
wards revealing dark matter's particle identity is to search for 
its signature in radiation backgrounds, either from the Milky 
Way or in the isotropic diffuse photon background (iDPB), 
which can contain both cosmological and Galactic halo con- 
tributions. Dark matter might be comprised of particles that 
can decay with finite lifetimes much longer than the age of 
the universe. In such scenarios, the resultant fluxes of decay 
products depend on the amount of dark matter present alone, 
as opposed to self-annihilation, which, being dependent on 
particle density squared, is very sensitive to assumptions con- 
cerning details of dark matter clustering. 

A wide variety of decaying dark matter models have been 
examined in regards to their observable implications 01- 
Among late decaying dark matter models, sterile neutrinos 
with multi-keV masses have been extensively studied as dark 
matter candidates PJ, with strong constraints placed on their 
decays, e. g. Refs. | H 0, Bjl and references therein. The de- 
cay of moduli dark matter f6j with masses of several hundred 
keV may contribute to the sub-MeV iDPB. The dark matter 
model of Ref. [7], inspired from minimal universal extra di- 
mensions or supersymmetry [8], with a mass scale of hun- 
dreds of GeV, is advocated as the source of the iDPB in the 
MeV range |9], which has yet to be accounted for with con- 
ventional sources (e.g., supernovae [lOJor active galactic nu- 
clei iHtl) or more exotic mechanisms I12I1 . Similarly, decay- 
ing gravitino dark matter in R-parity breaking vacua fl3], with 
multi-GeV masses, has been suggested as an explanation of 
iDPB spectral features in the GeV range il4ll5[l . 

Rather than focusing on a particular model, we first con- 
sider a generic decaying dark matter scenario in which the de- 
cay of the parent particle is dominated by a monochromatic 
photon emission. We assume that the lifetime of the par- 
ent particle, r, is much longer than the age of the universe 
(tq ~ 4.5 X 10^^ s), thus its cosmological abundance has 
not changed significantly since the time of dark matter decou- 
pling. The decay under consideration is x ~* x' + 7' where 7 
is a monochromatic photon emitted with energy e. In general. 



r and £ will depend on the masses of the parent and the daugh- 
ter particles (m^, m^') and their splitting. Am — — m^i. 

The flux of photons from dark matter decays is inversely 
proportional to both the particle lifetime (fixing the decay rate 
per particle as specified by a particular theoretical model) and 
the mass of an individual particle (yielding the total number of 
particles in a fixed amount of dark matter). Thus, gamma-ray 
observations allow us to place constraints only on the degen- 
erate product m^T versus e, as we display in Fig. [T] As we 
will discuss in detail, below the jagged line between 0.02- 
8 MeV, the gamma-ray line signal from the Galactic Cen- 
ter (GC) region due to dark matter decays violates the corre- 
sponding limit obtained with the SPI spectrometer on INTE- 
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FIG. 1: Model-independent constraints on the product of mass and 
lifetime, m^r, versus the energy carried away by the monochromatic 
photon emission, e, for a generic late-decaying dark matter model: 
X — > x' + 7- Regions excluded by either the gamma-ray line emis- 
sion limits from the Galactic Center region or overproduction of the 
isotropic diffuse photon background are shown, together with pre- 
ferred ranges of parameters from three well-studied models. 
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GRAL satellite ||16|]. Additionally, the iDPB, as determined 
from SPI ill, COMPTEL |18J and EGRET lUS ^ data, 
is overproduced (assuming it is fully accounted by late dark 
matter decays in a given energy band) in the triangular region, 
even disregarding any contributions from known astrophysi- 
cal sources. We also show three representative scenarios, in- 
spired by the theories of sterile neutrinos, R-parity breaking 
vacua, and mUED. Since m^, my and Am are not necessar- 
ily predetermined, they may be adjusted to yield the displayed 
curves relating m^r and e. 



n. MILKY WAY GAMMA-IfAY LINE SEARCH 

A monochromatic line will be most detectable locally, 
where cosmological redshifting is of no concern. Fortunately, 
a search for diffuse gamma-ray line emission in the energy 
range 0.02-8 Me V from the GC region has been conducted by 
Teegarden and Watanabe using the SPI spectrometer on the 
INTEGRAL satellite il6ll . which recovered the known astro- 
physical diffuse line fluxes, such as the 511 keV positron an- 
nihilation line 1I22I] . The excellent energy resolution of SPI en- 
abled them to place very strict constraints on potential uniden- 
tified emission lines, with an energy dependent 3.5 cr flux 
limit, !Fiim{E), from an angular region within a 13° radius 
of the GC (which we refer as the GC region). This limit, re- 
produced in the top panel of Fig. |2] can be compared to the 
expected gamma-ray flux arising from late dark matter decays 
in the GC region, which we calculate following the methods 
in Ref. L5J. 
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FIG. 2: Top: Limits on the diffuse gamma-ray line emission from 
the Galactic Center region (an angular region within a 13° radius) as 
adopted from Ref. [16]. Bottom: Representative measurements of 
the diffusephoton background from SPI 1 17], COMPTEL [ 18] and 
EGRET Hil] in the energy range around 0.01 MeV-100 GeV. The 
thick solid line, summarizing the overall trend of the data, is to be 
compared to predictions of decaying dark matter scenarios. 



We first define a dimensionless line-of-sight integral at an 
angle ijj relative to the GC, 



1 



PscRs 



2 e Rsc cos i; + e'^ 



(1) 

where p is the density of the dark matter in the halo as a func- 
tion of the distance from the GC. This is normalized to the 
dark matter density (p^c — 0.3 GeV cm^^) at the solar circle 
(Rg^ = 8.5 kpc) so that pscRsc ~ 8 x 10^^ GeV cm"^. Note 
that this arbitrary normalization is needed to make J' dimen- 
sionless and will be canceled-out later. The upper limit of this 
integration. 



2 

MW 



■ipR^.J + RscCosi/;, (2) 



depends on Rmw, the assumed size of the halo. J is rela- 
tively insensitive to £max as long as Rmw is large. The inten- 
sity of photons (number flux per solid angle) from the same 
direction. 



J(V.) = 



~Rs 



(3) 



can be integrated over a circle of radius around the GC (cov- 
ering a patch of area Aft = 2tt{1 — cos tp)) to obtain the cor- 
responding total flux. 
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dVL' Ji^,') (4) 



The limit reported in Ref. lfl6ll has been obtained by sub- 
tracting the average flux measured at regions away from the 
GC region (ip > 30°) from the average flux measured in- 
side the GC region (?/> < 13°) to eliminate instrumental back- 
grounds. Thus, the constraining power of this limit for decay- 
ing dark matter scenarios depends on the enhancement of the 
expected signal towards the GC region. Both theoretical and 
observational studies strongly suggest that the central regions 
of dark matter halos are significantly denser and, moreover, 
the column depth is higher towards the GC direction relative 
to off-axis lines-of-sight. We have reproduced the impact of 
this subtraction (see Ref. [5J for details) by calculating a pa- 
rameter 
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which ranges between ^ 0.5 — 1.5 for various dark mat- 
ter halo fitting profiles commonly used in the literature ll23ll . 
Here i7>3o° is the average of J' away from the GC region. 
We also note that the results that we adopted from Ref. ||T6ll 
are based on an assumption that the expected line signal has 
a Gaussian source profile, while a flat source profile could 
yield limits that are weaker by up to a factor of ~2 (e.g., see 
Fig. 12 of Ref. [|29|] ). IVIoreover, one would expect to see these 
limits improve as the amount of available data increases in 
time ll29ll . In the rest of our study, we choose a conservative 
value, ~ 0.5, which can be realized only for profiles that 
are rather flat inside the solar circle. While this mostly pro- 
tects our conclusions from uncertainties in the halo profile, our 
subsequent result can be easily rescaled for a different value. 
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The predicted gamma-ray emission line flux due to dark 
matter decays at a given e must not exceed the corresponding 
limits from the GC region, thus 



PscRsc 



Clim < J^\im{E = e) 



(6) 



Rearranging this equation yields our model independent con- 
straint. 



Psc^scCvim 



3 X 10^0 GeV cm" 

•?^lim(£) 



(7) 



as shown in Fig. [T] (region below the jagged line). The ex- 
pected dark matter decay flux is inversely proportional to 
m^T, which leads to an overproduction of gamma rays for 
iT^x^ ^ 10^^ ^^"^ s energy range 0.02-8 MeV. Thus 

the area below the jagged line is excluded by the the diffuse 
gamma-ray line emission limits from the GC region. 



III. ISOTROPIC DIFFUSE PHOTON BACKGROUND 



component is 



^ Psc-^SC rr 

-^iso — c/iso ■ 



(8) 



We present a representative spectrum for this isotropic sig- 
nal in Fig. [3] (dotted line), after convolution with a Gaus- 
sian of ^ 10% width to simulate the energy resolution of 
a typical detector We have chosen e — 1 MeV, with 
TO^r — 7 X lO'^"' GeV s, the maximum value allowed by the 
the line emission bounds from the GC region (Fig. [U. For 
these parameters, the isotropic contribution of the dark mat- 
ter decays in the Galactic halo alone to the iDPB is less than 
2% (in a bin of logarithmic width 0.4 dex centered around 
e — 1 MeV). Note that the average flux expected from the 
decays in the Galactic halo (which is more directional, peak- 
ing toward the GC region) can be at most several times larger 
than this isotropic component since we are dealing with de- 
caying dark matter particles (contrary to self-annihilating dark 
matter, which is highly sensitive to the details of dark matter 
clustering). 



While stringent limits on line emission from the GC re- 
gion are only available in a rather limited energy range (0.02- 
8 MeV), the iDPB is measured over a broad range by many 
instruments. In the bottom panel of Fig. |2] we display 
three recent determinations of iDPB in different ranges of 
energy from SPI d, COMPTEL iH and EGRET iH, 
which are consistent with others measurements (see, e.g., 
Ref. flli I24I1 ). The thick dotted line represents the global 
trend of the data to be used for comparison. We choose the 
terminology "isotropic diffuse" photon background (iDPB), 
as opposed to "cosmic" or "extragalactic", since the contribu- 
tion from sources in the Milky Way or its halo is not clear, and 
iDPB can include gamma-ray line signals that could not have 
been resolved by COMPTEL or EGRET. While it is gener- 
ally thought that AGN are responsible for the emission in the 

keV |25] and - GeV |26] ranges, the origin of the iDPB, 
especially in the MeV regime, is far from bein g se ttled, with 
various scenarios having been entertained [ fiol ll iL \l2\\ . It is 
then of interest to determine just how much of the iDPB can 
possibly be accounted for by late decaying dark matter 



A. Dark Matter Decays in the Halo 

While the photon signal from dark matter decays in the 
Galactic halo is enhanced towards the GC, as has been utilized 
for our constraints in the earlier section, it also contains an ap- 
parently isotropic contribution. The limited energy resolution 
of past gamma-ray detectors could not distinguish monochro- 
matic line emission from the Galactic halo from a truly cos- 
mological signal. The intensity of the isotropic halo contri- 
bution, 2iso, can be estimated from a line of sight integration 
in the anti-GC direction, Ji^o = ^7(180°) 1, as this is the 
minimum contribution from the dark matter halo of the Milky 
Way. Regardless of the underlying halo profile, this number 
is relatively robust, being mostly dependent on the dark mat- 
ter density at the solar circle. The intensity of this isotropic 



B. Cosmological Dark Matter Decays 

We now evaluate the contribution of truly cosmological 
dark matter decays to the iDPB. For late decaying particles 
(t ;2> To), the comoving dark matter density has remained 
nearly constant since the early universe. The comoving de- 
cay rate is then simply proportional to the dark matter fraction 
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FIG. 3: Photon spectrum from isotropic Galactic halo decays (dotted 
line) for e = 1 MeV, with m^T ~ 7 x 10^* GeV s chosen from 
Fig. [T] such that the line emission bounds from the Galactic Center 
region are saturated. Also displayed are the spectra from cosmolog- 
ical decays (dashed line) and the total spectrum (solid line), which 
falls well below the isotropic diffuse photon background (thick solid 
line). 
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~ 0.25) of the critical density of the universe, pc, and 
is given as pc^^/ {ra^T). The diffuse gamma-ray flux (per 
solid angle per unit energy) arising from the decays can be 
calculated by considering the contributions from all redshifts 
(analogous to L27ll ). 

ah An Hq J n{z) rriy^r 

where h{z) = [(1 + z^Qm + ^a]^^^, ^ 0.3, Qa ~ 0.7, 
Ho = 70 kms"^Mpc~\ and c = 3 x IQi'' ems^^ (so that 
c/Ho ~ 1.3 X lO^Scm and p^ = 5.3 x lO'^ GeV cm^^). 
The integration can be eliminated after using the (5-function 
identity; S{ax — 6) = 5{x — b/a)/a, simplifying the result to 

d$ _ 1 e p^n^ 1 e(£ - E) 

dE 47r Ho m^T E ^{e/EY^Im + f^A ' 

where h is substituted and is a step function. We show this, 
using the same parameters as in the preceding subsection and 
again smoothing with a ~ 10% Gaussian, in Fig. [3] (dashed 
line). As seen in the figure, this cosmological flux is slightly 
lower than the isotropic contribution from the Galactic halo, 
and their sum (solid line) still falls well short of the observed 
signal, restricting their combined contribution to the iDPB to 
be less than 4%. 

To quantify and generalize our observations, we calculate 
the expected total (cosmological plus isotropic Galactic halo) 
spectrum for all values of m^r and compare to the iDPB (as 
denoted by the thick trend curve in Fig. |2]l, integrating both 
in a bin of logarithmic width 0.4 dex centered around e. This 
choice encompasses most of the expected signal where the 
decay spectrum peaks, and both exceeds the experimental en- 
ergy resolution and the uncertainties on the determination of 
the iDPB. In Fig. [T] the region in which dark matter decays 
overproduce the iDPB is shown (triangular region). 

Above this region, decaying dark matter alone cannot fully 
account for the iDPB. In fact, since there should be addi- 
tional contributions from AGN at both low and high ener- 
gies L25> .26.1 , the actual bound on the parameter m^r will 
be even more stringent than the one presented. Combining 
the iDPB overproduction constraint and the gamma-ray line 
emission limit from the GC region model-independently ex- 
cludes a sizable region in the parameter space of m-^T versus 
e, with the latter picking up when the former is exhausted at 
e ~ 8 MeV. 



IV. DECAYING DARK MATTER MODELS 

While we derive our constraints for a decay scenario that 
is dominated by monochromatic photon emission, there may 
be additional modes of decay or self-annihilations producing 
other signals. Our constraints on the lifetime of the dark mat- 
ter candidate via monochromatic photon emission could be 
generalized to the total lifetime including other decay chan- 
nels, as long as the latter is long enough to justify the assump- 
tion that the cosmological abundance of the parent particle has 
not changed significantly. 



For the generic decay we are considering, the energy of the 
emitted photon is dictated by the splitting. Am, as follows. 
When Am <C m^/ (or equivalently ~ n^x')' recoil 
of the daughter can be neglected, so that e — + Am. For 
Am ^ m^i (or m^ ^ ni^'), two relativistic particles are 
produced, so that e Am/2 ~ "^x/-^- Generally, models lie 
in one of these two regimes. To emphasize the generality of 
our constraints, now we discuss particular scenarios. 

For example, WlMPs with weak-scale masses and cross 
sections may have monochromatic decays. The decay process 
between the two lightest Kaluza-Klein (KK) particles, the KK 
hypercharge gauge boson, and KK graviton, in mUED 
models, and the decay between the two lightest particles in 
SUSY theories, the Bino-like neutralino, B and gravitino, G 
are well-studied examples. The mass scales of these candi- 
dates are --800 GeV for the former and --80 GeV for the lat- 
ter. The decay rates in these theories [9] are highly suppressed 
due to the weakness of gravity and is given by 

4.7 X 10^2 s / Am \ 

b (m^J ■ ^''^ 

where the parameter b is identified as (2, 10/3, 1, 2) for each of 
the decay reactions (G^ B^+j, B^ G B+7, 

B ^ G + j). The lifetime requirement of t ^ tq translates 
into Am < 30 MeV. Since Am ^ m^, the energy carried 
away by the emitted photon is e ~ Am. Eq. (fTTI) can be 
rearranged as 

™,..4.7xl0-.(!ii)(^)-\ ,12) 

which relates m^r to e in terms of a single parameter: m^/b. 
We plot m^T versus e in Fig.|4]for m^/b ~ 300 GeV to repre- 
sent mUED. One sees that the Milky Way constraint requires 
e < 1.5 MeV, which is a very strict limit as the lifetime is 
proportional to e^'', i.e., the decay rate increases by almost an 
order of magnitude from 1 MeV to 2 MeV. This translates 
to the restriction of Am < 1.5 MeV, which is far stricter 
than the necessary condition to have a long-lived candidate. 
Am < 30 MeV. 

In the R-parity violating supersymmetric extension of the 
standard model, the lightest supersymmetric particle is again a 
gravitino that might not be stable on cosmological time scales 
against decay into a photon and neutrino (G ^ 1^ + 7) through 
a small photino-neutrino mixing \U^,y\. The lifetime of the 
gravitino in this model ifTsIl is 

with the resulting photon and neutrino each carrying an energy 
of e = '7ijj,/2. We can rewrite this equation in terms of m^r 
versus e as 

mxT :^ 1014 GeVs (|C/^.r' (^)"' . (14) 

We plot this relation for \U^,y\ = 10"* in Fig. [T] which 
shows that the contribution of this decay model to the iDPB 
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FIG. 4: Similar to Fig. [T] focusing upon the MeV range of e. The 
contribution of late dark matter decays to the isotropic diffuse photon 
background is 10% or more in the diagonal band. SUSY and mUED 
inspired decaying dark matter models of Ref . [9J cannot make signif- 
icant contribution to the iDPB while abiding by the ganraia-ray line 
emission limits from the Galactic Center region. 



will be significant around e ~ 5 GeV (corresponding to 
~ 10 GeV) agreeing with Ref. ifTsll . Slightly above/below 
^ IQ GeV, either its contribution is negligible or vastly 
overproduces the iDPB. 

Dark matter models involving keV-mass sterile neutrinos, 
in their simplest description, require only two parameters, the 
sterile neutrino's mass and mixing with active neutrinos. The 
decay chain for sterile neutrinos is Vg Vc,fi.T + 7, with a 
radiative lifetime Ii28il (for Dirac neutrinos) of 

T = 1.5 X 10^2 s(sin"2 261) \ (15) 

This can similarly be rearranged, keeping in mind that the en- 
ergy of the parent sterile neutrino is split equally between the 
photon and the daughter neutrino (e = ms/2), as 

rUsT ~ 10^5 GeVs (sin^^ 261) , (16) 

which has only one free parameter, sin'^ 26*. For illustration, 
we plot Eq. ( fT^ in Fig.[T]for sin^ 26 = 10~^^. As seen in the 
figure, and has been established in Ref. ifsl l29ll . the gamma- 
ray line emission limit from the Galactic Center region pro- 
vides quite stringent restrictions on sterile neutrino dark mat- 
ter, which can be several orders of magnitude stronger than 
constraints from overproduction of the iDPB. Interestingly, all 
three models we have discussed have the form m^r oc e~", 
where a = 3, 2, 4 respectively, as can also be noticed from the 
varying slopes of the lines representing the models in Fig.II] 



V. CONCLUSIONS 



Predictions of photon fluxes from dark matter decays are 
considerably more robust than those from annihilations, due 
to a lesser dependence upon theoretical uncertainties in the 
distribution and clustering of dark matter We have shown 
that the gamma-ray line emission limits from the Galactic 
Center region, along the isotropic diffuse photon background, 
allow for stringent constraints to be placed on late decaying 
dark matter scenarios that produce monoenergetic photons. 
We emphasize that the Galactic and cosmic constraints are 
not independent of each other, with the GC region provid- 
ing stronger limits in its range of applicability due to new 
spectroscopic data. Rather than attempting to explain vari- 
ous gamma-ray phenomena with a specific model, we report 
model-independent bounds on the decaying dark matter pa- 
rameter space (as defined by m^r versus e). Our general 
constraints are applicable to a number of models, and can 
be used as a guide for future model building. Upcoming 
gamma-ray telescopes with improved energy resolution, such 
as GLAST |30] or ACT IsTll . can improve upon these bounds, 
particularly by making use of the unique spectral shape and 
directionality of decays from the Galactic halo. 

One interesting application of our study is to assess the 
recent suggestion that cosmological late dark matter decays 
can explain the isotropic diffuse photon background in the 
1-5 MeV range, whose origin remains a mystery f^. We 
plot m^T versus e in Fig. H] for m^/h ~ 300 GeV and 
m^/b ~ 50 GeV to represent the aforementioned mUED and 
SUSY models of Ref. i^, respectively. We also show the 
range of parameters, m^r versus e, that can lead to a substan- 
tial (> 10%) contribution to the iDPB (shaded diagonal band) 
or overproduce them (triangular region) through the sum of 
the local decays (Galactic halo) or decays from truly cosmo- 
logical sources (all distant dark matter halos). The region ex- 
cluded by the gamma-ray line emission limits from the GC 
region is below the jagged line. As seen in the figure, even the 
combined emission from the Galactic halo and cosmological 
sources due to either the mUED or SUSY inspired decaying 
dark matter models of Ref. |9] cannot make a significant con- 
tribution to the iDPB while abiding by the gamma-ray line 
emission limits from the GC region. The mUED model can 
contribute to the iDPB only for Am < 1.5 MeV with a contri- 
bution of < 5%, while the SUSY model is even more severly 
constrained. Even relaxing our assumptions on the distribu- 
tion of dark matter in the halo does not increase these frac- 
tions dramatically, thus, dark matter cannot decay in the late 
universe at a high enough rate to make a prominent contribu- 
tion to the iDPB in the MeV range. 
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